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ABSTRACT 


In  this  report,  fatigue  in  general  and  some  prominent  theories 
concerning  prediction  of  cumulative  fatigue  damage  are  discussed.  A 
computer  program  was  developed  to  calculate  the  cumulative  fatigue 
damage  and  fatigue  life  using  the  predictive  equation  developed  by 
I.  R.  Kramer  (8).  Test  results  generated  by  Kramer  for  2014-T6 
aluminum  alloy  were  used  to  determine  cumulative  fatigue  damage  and 
fatigue  life.  The  experimental  values  of  fatigue  damage  and  life  are 
found  to  be  in  agreement  with  those  predicted.^ 
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PREDICTION  OF  CUMULATIVE  FATIGUE  DAMAGE 


Failure  of  machine  members  under  repeated  or  fluctuating 
stresses  is  called  fatigue  failure.  This  type  of  failure  occurs 
below  the  ultimate  strength  of  the  material  and  quite  often  even  below 
the  yield  strength.  Failure  begins  with  a  crack  at  the  surface.  The 
initial  crack  is  so  small  and  minute  that  it  cannot  be  detected  by  the 
naked  eye  and  is  even  diffivcult  to  locate  in  X-ray  inspection.  The 
crack  developes  at  a  point  of  discontinuity  such  as  a  keyway,  a  hole, 
an  inspection  or  stamp  mark,  an  internal  crack,  or  some  irregularities 
caused  by  machining.  Once  a  crack  has  developed,  the  stress  concentra¬ 
tion  effect  becomes  greater  and  the  crack  propagates  more  rapidly.  As 
the  stress  area  decreases  in  size,  the  stress  increases  in  magnitude 
until,  finally,  the  remaining  area  fails  suddenly.  A  fatigue  failure 
is,  therefore,  characterized  by  two  distinct  areas  of  failure.  The 
first  is  due  to  the  progressive  development  of  crack,  while  the  second 
is  due  to  sudden  failure  which  resembles  the  fracture  of  brittle  materials. 
Many  static  failures  are  visible  and  give  warning  in  advance.  But  a 
fatigue  failure  gives  no  warning,  it  is  sudden  and  total,  and  hence 
dangerous. 


Fatigue  Strength  and  Endurance  Limit 


To  establish  the  fatigue  strength  of  a  material,  quite  a  number 
of  tests  are  necessary  because  of  the  statistical  nature  of  fatigue.  The 
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first  test  is  made  at  a  stress  which  is  somewhat  under  the  ultimate 
strength  of  the  material.  The  second  test  is  made  with  a  stress  which 
is  less  than  that  used  in  the  first.  This  process  is  continued  and  the 
results  are  plotted  as  an  S-N  diagram  (Fig.  1).  This  chart  may  be 
plotted  on  semi-log  or  on  log-log  paper.  In  the  case  of  ferrous  metals 
and  alloys,  the  graph  becomes  horizontal  after  the  material  has  been 
stressed  for  a  certain  number  of  cycles.  The  ordinate  of  the  S-N 
diagram  is  called  the  Fatigue  Strength,  corresponding  to  the  number  of 
cycles  N  required  to  produce  failure.  When  the  curve  becomes  horizontal, 
as  it  does  for  steel,  failure  will  not  occur  if  the  stress  is  below  this 
level,  no  matter  how  many  stress  cycles  are  applied.  This  fatigue 
strength  is  called  the  Fatigue  Limit  or  Endurance  Limit.  Different 
stress  components  used  in  fatigue  analysis  are  shown  in  Figure  2. 

Low-Cycle  and  High-Cycle  Fatigue 

A  complete  S-N  curve  may  be  divided  into  two  portions;  the  low- 

cycle  range  and  the  high-cycle  range.  There  is  no  dividing  line  between 
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the  two.  The  investigators  arbitrarily  say  that  up  to  about  10  or  10 

4 

is  low  cycle  and  beyond  10  cycles  is  high  cycle.  The  low-cycle  fatigue 
is  importnat  in  pressurized  fuselages,  missiles,  space  ship  launching 
equipment,  etc.  The  failure  mechanism  in  the  low-cycle  range  is  close 
to  that  in  static  loading,  but  the  failure  mechanism  in  the  high-cycle 
range  is  different  and  may  be  termed  "true  fatigue." 

Cumulative  Fatigue  Damage 

Fatigue  loads  applied  to  machine  parts  and  structures  are 
seldom  of  constant  magnitude.  Machines  have  to  be  started  up  and 
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stopped,  overloads  occur,  and  transient  vibrations  of  the  part  may 
impose  high  frequency  stresses. 

In  Figure  3  assume  that  a  specimen  is  subjected  to  an  alternating 
stress  of  amplitude  for  cycles.  The  fatigue  strength  corre¬ 
sponds  to  Ng  cycles  of  life,  and  the  remaining  life  at  this  same  stress 
magnitude  is  Ng  -  cycles.  Consequently,  the  specimen  has  accumulated 
some  fatigue  damage  at  this  stress  magnitude.  OBD  represent  the  S-N 
diagram  of  the  virgin  specimen.  =  Ng  -  is  the  remaining  useful 
life  of  the  specimen.  Now  locate  point  C,  and  construct  line  OCE,  which 
is  the  new  S-N  diagram  having  a  lower  endurance  limit  o'.  The  damage 
done  by  overstressing  is,  therefore,  the  difference  in  the  endurance 

limits  (a  -  a'). 
e  e 

Equivalent  Cycle  Approach.  Let  us  consider  the  two-level 
sinusoidal  stress  history  shown  in  Figure  4  with  maximum  stress  levels 
S-j  and  $2  applied  to  a  material  alternately  in  groups  of  n"^  and  172 
cycles,  respectively.  The  number  of  cycles  to  failure  at  stress  condi¬ 
tions  S-j  and  $2  are  represented  by  N-j  and  N2  respectively.  From 
Figure  5; 


or 

-  Ik- 

"21  ■  "1 

n2i  is  the  number  of  cycles  applied  at  stress  condition  $2  that  would 
produce  the  same  amount  of  damage  as  77^  cycles  applied  at  stress  condition  S^. 


By  a  similar  analysis 


^12  "2 

The  number  of  cycles  n2g  at  stress  condition  that  yields  the 
same  amount  of  damage  as  that  caused  by  the  block  containing  n^ ,  cycles 
at  stress  condition  S-j ,  plus  n^  cycles  at  stress  condition  $2- 

^2B  "  '"'21  ^  ^2 


=  N 


2^N, 


The  number  of  repetitive  blocks  to  failure  n^^  is 


n 


Bf 


"2B 


n 


Bf 


Replace  n2g 


n,  n  n  n 


or 


2  n  , 
y  =  1 
S  N 
q=l  q 


Equivalent  Stress  Approach.  The  equivalent  stress  approach 
derives  its  name  from  the  consideration  that  there  exists  a  stress 
condition  which  will  cause  failure  in  the  same  total  number  of  cycles 
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as  that  needed  by  the  complex  history.  Total  damage  at  failure  Dp, 
associated  with  a  multilevel  h  distinct  sinusoidal  group  is  given  by 


V 

N 

q=l  q 


where  n  r  is  the  number  of  cycles  in  the  failure  history  at  stress  condi- 

qF 

tion  S  .  The  total  number  of  cycles  required  to  cause  failure  is  Np. 

q 

Equivalent  stress  approach  specifies  that  a  stress  condition  would 
also  produce  failure  in  Np  cycles.  If  Dp  denotes  the  damage  at  failure 

associated  with  the  stress  condition  S  then  D^.  =  D^  . 

^  ^  '’e 

If  Ng  is  the  number  of  cycles  to  failure  at  stress  condition 

then 


N 

e  e 


THEORIES  ABOUT  THE  PREDICTION  OF 
CUMULATIVE  FATIGUE  DAMAGE 


Miner's  Theory 

The  simplest  and  most  often  used  is  the  theory  proposed  by  Miner. 
This  theory  is  referred  to  as  the  linear  cumulative  damage  rule  and 
utilizes  the  simple  cycle  ratio  as  its  basic  measure  of  damage.  If  a 
multilevel  sinusoidal  stress  history  is  applied  to  a  structural  material, 
it  is  hypothesized  in  Miner's  Theory  that: 

a)  each  group  of  sinusoids  contributes  an  amount  of  damage 
given  by  the  linear  cycle  ratio  for  the  group. 

b)  the  damage  done  by  any  group  of  sinusoids  is  not  dependent 
on  the  group's  location  in  the  stress  history. 

c)  the  total  applied  damage  is  equal  to  the  sum  of  the  damages 
contributed  by  each  sinusoidal  group. 

n 

D  =  J- 
q 

D^  =  The  damage  resulting  from  this  group  of  sinusoids, 
n^  =  Number  of  sinusoids  in  this  group,  and 
Nq  =  Number  of  sinusoids  to  produce  failure  at  maximum  stress 
level  S  . 

q 

If  Dg  is  damage  produced  by  the  block  of  h  distinct  sinusoidal 
groups  then 


6 
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D 


B 


h 

I 

q=l 


n 

_9. 

N 


If  there  are  n„  number  of  the  basic  block  of  h  sinusoidal  group, 

D 

the  total  damage  D  is  given  by 


h 


n  h 

N  ^ 

q  q=i 


n 

N 


For  failure  to  occur  D  must  equal  unity. 

It  has  been  found  from  testing  under  multilevel  sinusoidal 
histories  that  Miner’s  theory  predicts  a  longer  life  than  that  actually 
wi tnessed. 


Grover's  Theory 

Grover's  theory  considers  that  the  fatigue  life  of  a  material 
subjected  to  a  complex  stress  history  is  composed  of  two  stages:  a)  an 
initial  number  of  cycles,  Np  required  to  nucleate  the  crack  and  b)  iNp  num¬ 
ber  of  cycles  needed  to  propagate  this  crack  to  the  failure  of  the  material. 
The  total  number  of  cycles  Np  required  to  cause  the  failure  of  the  material 
is  thus  given  by 

Np  =  N’  +  Njl 

Now  we  consider  a  multilevel  sinusoidal  failure  hissory  with  h  distinct 
sinusoidal  stress  conditions 


N 


F 


and 


Njl 


h 

I 

q=l 


n 


qF 
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where  n^p  and  n^p  designate  respectively  the  number  of  cycles  at  stress 

condition  S  applied  during  the  crack  nucleation  and  crack  propagation. 
Q 

If  is  the  number  of  cycles  in  the  failure  producing  crack  and  is 

the  number  of  cycles  needed  to  propagate  the  crack  to  failure  at  stress 

S  and  N  is  the  number  of  cycles  to  failure  at  S  then 
q  q  q 

N  =  N'  +  N" 

q  q  q 


N'  N" 

+  _a  = 

N  N 

q  q 


Grovers  theory  now  utilizes  Miner’s  theory  separately  for  the 
nucleation  stage  and  for  the  progation  to  failure  stage. 


y  = 


y  ^  = 
Z  N" 


This  theory  has  a  serious  setback  of  inability  to  find  exactly 
the  number  of  cycles  required  to  nucleate  the  crack.  This  theory  is 
unconservative  like  Miner's  Theory. 


Marco-Starkey  Theor\ 


The  Marco-Starkey  specification  for  the  damage  0  arising  from  n 
cycles  applied  at  stress  condition  S  with  an  associated  number  of  cycles 
to  failure  N  is  given  by: 

D  =  (n/H)^ 
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The  exponent  x  is  a  variable  quantity  whose  magnitude  is  dependent 
on  the  applied  stress  condition.  Marco  and  Starkey  consider  that  x  has 
a  magnitude  greater  than  unity  and  approaches  unity  as  the  stress  condi¬ 
tion  becomes  severe  as  shown  in  Figure  6.  This  theory  is  conservative 
and  has  limited  use  because  of  the  difficulty  to  determine  the  exponent  x, 
and  its  dependency  on  stress  under  complex  cyclic  conditions. 


n  =  number  of  cycles  applied  at  stress  condition  S. 
b  =  the  slope  of  the  central  portion  of  the  S-N  diagram. 

C  &  k  =  material  constants  where  k  is  greater  than  1. 

D  is  a  function  of  number  of  applied  cycles  rather  than  the  cyclic 
ratio  n/N  as  in  previous  cases. 

It  is  seen  that  the  equation  for  the  central  portion  of  the  S-N  diagram 
can  be  put  in  the  form 

N  =  , 

C 

therefore , 


N  and  S  are  reference  number  of  cycles  and  stress  respectively.  If 
r  r 

we  take  k  =  1  then  D  =  n/N,  which  is  Miner's  theory.  The  value  of  k  =  2 
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is  mostly  used  for  Shanley's  theory.  This  theory  yields  shorter  fatigue 
life  than  that  predicted  by  Miner. 

Corten-Dolan  Theory 

For  the  Corten-Dolan  theory  the  damage  0  in  a  material  due  to 
n  cycles  of  pure  sinusoidal  stress  history  can  be  expressed  as: 

D  =  r  n^  (1) 


r  =  function  of  stress  condition 
a  =  material  constant 

kb 

If  a  =  1  and  r  =  C  S  then  Corten-Dolan  theory  yields  to 
Shanley's  theory.  Corten  and  Dolan  used  an  equivalent  cycle  approach 
for  determining  a  fatigue  failure  criterion  based  on  the  damage  specifi¬ 
cation  in  Eq.  ( 1 ) . 


Freudenthal-Hel ler  Theory 

Freudenthal  and  Heller  proposed  the  modification  to  the  central 

portion  of  the  S-N  diagram.  In  terms  of  modified  number  of  cycles  to 

failure  N  for  stress  condition  S  can  be  expressed  as 
qm  q 


N  S* 

-511  =  (^)6 
N*  ’ 

r  q 


(2) 


where  S*  =  reference  stress  condition  that  is  unrelated  to  any  applied 
stress  history, 

N*  =  associated  number  of  cycles  to  failure,  and 
6  =  slope  of  modified  S-N  diagram. 


n 


Expression  for  the  conventional  S-N  diagram 


_9.  = 

N* 


(3) 


The  above  equations  (2  &  3)  indicate  that  the  modified  and  conventional 

S-N  diagrams  will  be  coincident  only  at  stress  condition  S  =  S*  for 

q  r 

which  Nq  =  N*.  The  reference  stress  condition  S*  will  usually  have  such 
intensity  that  N*  will  fall  in  the  range  of  10^  to  10^  cycles.  Total 
number  of  cycles  N  required  to  cause  failure  of  a  material  under  a 

multilevel  sinusoidal  stress  history  is  as  follows,  based  on  the  modified 
S-N  diagram  shown  in  Figure  7. 


N 


h 

F  ^ 

m  q=l  qm 


The  total  number  of  cycles  Np  required  to  cause  failure,  based  on 
the  conventional  S-N  diagram,  is  expressed  as 


rh  1 

rh 

a 

- j  n 

_ 

iq=l 

i_ 

_9_ 

N 

h 

V 

I 

q=i 

h 

y 

L, 

q=l 


Kramer's  Theory 

Kramer  conducted  some  experiments  and  proposed  that  while 
materials  are  subjected  to  fatigue  cycles,  the  work  hardening  of 
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surface  layer  takes  place  and  consequently  the  proportional  limit  for 
the  material  is  increased  with  increased  number  of  cycles.  He  defined 
this  increase  in  proportional  limit  as  the  surface  layer  strength  (a^) 
He  further  investigated  that  when  this  surface  layer  stress  reaches  a 
critical  value  (o*)  the  failure  producing  crack  is  propagated.  He 
showed  that  is  independent  of  the  stress  magnitude  applied. 

Nq  =  number  of  cycles  to  initiate  the  propagating  crack 
Np  =  number  of  cycles  to  failure 
N 

ir-  =  Const  =0.7  for  aluminum 

"f 

S  =  d^/dN 

a  =  SN 
s 

or  a*  =  S 
s  0 

a  ,  °s. 

D  =  -T  a  crack  will  initiate  when  T  =  1  or  S.N.  =  a* 

O*  “  IIS 


The  incremental  rate  of  change  of  surface  stress  at  the  first 
stress  level  is  given  by: 

S,  =  a  P, 

1  o  1 


After  cycles,  the  maximum  stress  is  increased  to  and  the 
incremental  rate  of  change  of  surface  stress  at  this  second  level  will 
be  modified  as 


S 


II 


f 

J 
S 


P 

^2 


J 


Similarly  at  third  stress  ievel 


■ill 


(57*  “3 


Replace  S 


II 


III 


0,  "1  4  '2  p 

'o  )  *  p)  "3 

^  “3 


III  '0, 


p,  'fifa  o,  "2  p 
(5r>  <>3 


And  so  on. 


S,N, 


=  n* 


p  p  a,",  .  c/f,f2p/ 

“i"!  *  *  “3"3(7>  (7) 


“s'k,  o,  “3^N,  p,  "1^2 

*  V4> 


^1  ’  ^2’  of  applied  cycles 

Ol .  O2,  o^'-applied  stresses 
^1  ’  ^2’  ^2“"P'"6''ious  history  damage  terms 
P  =  -1/m  — m  is  slope  of  S-N  curve 
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C  =  material  constant 

log  0^  -  log  0^ 

"  log  -  log 

Log  o  =  mLog  N  +  Log  C 
=>  C 

Log  a  =  Log  +  Log  C 

Log  0  =  Log  C  N*” 
or 

a  =  C  N*"  =c.  N*"  =  I 

N  »  o’^”  .  C-’-'” 

N  •  c-f  C" 
or 

-P 

N  =  e  0 

Kramer's  equation  uses  only  the  S-N  diagram  and  also  takes  care 
of  the  previous  damage  histories  in  the  following  terms. 

Experimental  Work.  The  fatigue  specimens  used  to  measure  the 
surface  layer  stress  and  to  determine  the  effects  of  removing  the  surface 
layer  on  fatigue  life  were  machined  from  15mm  diameter  rods  of  2014-T6 
aluminum.  These  specimens  had  a  diameter  of  0.16  and  a  gage  length  of 
0.30  in.  Before  testing,  all  specimens  were  electropol ished  to  remove 
about  0,004  in  to  obtain  the  same  surface  finish.  The  fatigue  tests 
were  conducted  in  an  electrohydraulic  machine  in  tension-compression. 
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The  change  in  the  surface  layer  stress  at  various  applied  stress 
amplitudes  was  determined  by  measuring  the  increase  in  the  proportional 
limit  as  a  function  of  the  number  of  cycles.  However,  when  the  surface 
layer  was  removed  after  cycling,  the  proportional  limit  decreased  to  the 
same  value  as  that  of  the  unclycled  specimen.  It  follows  from  these 
fatigue  tests  that  the  work  hardening  of  the  specimens  during  cycling  is 
confined  primarily  to  the  surface  layer.  The  increase  in  the  proportional 
limit  is  then  equal  to  the  strength  of  the  surface  layer. 

For  the  determination  of  the  change  in  the  proportional  limit, 
the  stress-strain  measurements  were  made  immediately  after  cycling  to 
minimize  surface  layer  losses  due  to  relaxation  effects.  An  extensometer 
with  gage  length  of  0.30  in  was  attached  to  the  specimen  after  the 
cycling  sequence  and  measurements  were  begun  in  less  than  one  minute. 

Analysis.  To  verify  the  validity  of  his  equation,  Kramer  fatigued 
seven  specimens,  using  four  stages  of  stress  with  different  combinations 
of  number  of  cycles  employed.  He  used  four  different  stress  sequences, 
i.e.,  low  to  high,  low  to  high-mixed,  high  to  low,  and  high  to  low-mixed. 

In  Table  1,  specimen  numbers  1  through  3  are  stressed,  using 
low  to  high  stress  pattern,  numbers  4and5  are  stressed  high  to  low, 
number  6  is  tested  low  to  high-mixed  and  number  7  is  tested  following 
high  to  low-mixed  stress  sequence.  Table  2  compares  the  total  cumulative 
fatigue  damage  calculated  using  Kramer's  and  Miner's  equations.  Table  3 
shows  actual  number  of  cycles  taken  by  specimens  in  the  final  (3rd  or  4th) 
stage  and  the  number  of  cycles  predicted  using  Kramer's  equation  in  the 
final  stage. 


CONCLUSIONS 


1.  The  predicted  life  is  in  good  agreement  with  that  determined 
experimentally  for  reversed  stress  conditions. 

2.  The  material  constants  P  and  0  are  calculated  from  S-N  diagram. 

The  accuracy  in  predicting  the  fatigue  life  under  cumulative  damage 
depends  to  a  large  degree  on  P  and  6.  P  and  3  vary  tremendously 
with  the  slope  of  S-N  diagram.  Therefore,  the  accurate  generation 
of  S-N  diagram  is  the  most  important  factor  for  obtaining  good 
results  using  Kramer's  equation. 

3.  The  cumulative  damage  values  calculated  by  Kramer  and  computed  at 
Tuskegee  Institute  differ  a  little  bit.  Most  probably  explanation 
could  be  that  Kramer  rounded  off  the  experimental  stress  values. 
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Figure  2.  Fatigue  Stress  Components 
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TABLE  1 

ACTUAL  EXPERIMENTAL  STRESS  SEQUENCE 


Specimen 

No . 

Stress 

Stress 

in  MPA/Nos.  o 

Sequence 

f  Cycles  in  K-Cycles 

1 

133/50 

172/40 

207/25 

241/9.2 

(F) 

0 

138/100 

172/75 

207/24 

(F) 

3 

133/25 

172/30 

207/24 

276/6.6 

(F) 

4 

276/5 

241/10 

207/24 

(F) 

5 

276/10 

241/5 

207/20 

172/56 

(F) 

6 

138/50 

241/10 

172/30 

276/10 

(F) 

7 

241/10 

172/30 

276/5 

207/27 

CF) 

(F)  Indicates  failure. 


TABLE  2 


ICIMEN 

1 

1 

CUMULATIVE 

DAMAGE 

KRAMER 

:  COMPUTED 

1 

1.1 

1.2961 

2 

1.4 

1.5267 

3 

1 

1.0 

1.0623 

1.1 

1.0744 

5 

1.8 

1.8779 

6 

1.1 

1.1550 

7 

1.4 

1.4366 

TABLE  3 


NUMBER  OF  CYCLES  TO  FAILURE 


SEPCLMEN 


EXPERLMENTAL 

PREDICTED 

1 

9200 

9000 

2 

24000 

29000 

3 

6600 

7000 

24000 

27  200 

5 

56000 

22200 

6 

10000 

20000 

27000 


39000 
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